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Abstract 

Heavy Quark Effective Theory splits a heavy quark momentum 
into a large fixed momentum and a variable residual momentum, = 
niQVfj, + kfj.. It thereby suffers a redundancy of description correspond- 
ing to small changes in the choice of the fixed velocity, u^. The fact 
that full QCD is manifestly ^-independent should lead to a non-trivial 
constraint on the form of the effective theory, known as Reparame- 
terization Invariance. For spin-1/2 quarks, the precise form of the 
constraint and its solution at the level of the effective lagrangian have 
proven to be rather subtle, and the original proposal by Luke and 
Manohar has been questioned. In this paper I employ a version of 
Heavy Quark Effective Theory containing the "anti-particle" field as a 
non-propagating auxiliary field, which greatly simplifies keeping track 
of i;-dependence. This permits a very simple derivation of Reparame- 
terization Invariance from first principles. The auxiliary field can also 
be integrated out to return to the standard formulation of the effective 
theory, but with the effective lagrangian now satisfying the full repa- 
rameterization constraint. I compare this result with earlier proposals. 
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1 Introduction 



Heavy Quark Effective Tlieory (HQET) provides a systematic expansion of 
QCD amplitudes in powers of (momentum transfer)/m, where m is a heavy 
quark mass. Calculations are very conveniently done in terms of an effective 
lagrangian which is also organized in powers of \/m. Ref. ||l| contains a 
review and references to the original literature. There are many interesting 
processes where the momentum transfers between a heavy quark and the 
light quarks and gluons are of order Aqcd- For such situations and for 
m S> AqcDi the heavy quark velocity is approximately constant. This is 
formally reflected in setting up the 1/m-expansion by splitting the heavy 
quark momentum into two pieces, 

where is a fixed four- velocity {v^ = 1), and A;^ is a variable residual 
momentum of order Kqcd- While this split is very useful for book-keeping, 
clearly there is an arbitrariness of order KqcD/'m in the precise choice of 
for any process [||. Formally we need just note the redundancy of the 
momentum decomposition under infinitesmal shifts in the fixed velocity, 

v' = V + 5v, v.5v = 0. (2) 

QCD is manifestly invariant under such infinitesimal reparameterizations, 
so a systematic approximation to full QCD such as HQET must also re- 
flect this fact. This constraint on HQET is known as Reparameterization 
Invariance (RPI). Luke and Manohar Q proposed that RPI should directly 
constrain the form of the lagrangian, like any other symmetry. Since the 
HQET lagrangian must be determined order by order in ag and 1/m by 
matching to full QCD, there is a significant reduction of work and increase 
in understanding if one can first enforce RPI on the general form of the 
effective lagrangian. 

The precise form of the RPI constraint on the HQET lagrangian for 
spin-1/2 quarks has been rather difficult to identify and prove. The origi- 
nal proposal by Luke and Manohar Q for the form of the RPI-constrained 
lagrangian was shown by Chen to be inconsistent with the direct match- 
ing computation already at tree-level. It is possible that the proposal is 
valid up to field redefinitions (see ref. Q for an application of RPI by this 
means) but this is less useful. Chen proposed another form for the RPI 
constraint consistent with tree-level matching, but its general validity has 
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been unclear because the derivation effectively depends on the use of the 
classical equations of motion. Kilian and Ohl obtained a RPI form for 
the HQET lagrangian valid to all orders, and which reduces at tree- level to 
Chen's version. I believe this to be correct and part of my present purpose 
is to give a simpler derivation and better understanding of this result. RPI 
in HQET has most recently been investigated by Finkemeier, Georgi and 
Mclrvin Q, who showed that the constraints from Chen's version of RPI 
are correct to order 1/m?. The framework of the present paper will clarify 
why this happens even though Chen's RPI is wrong in general (see section 
6). Ref. also showed that there is a field redefinition (which does not 
affect the S-matrix) taking Chen's constraints to those predicted by Luke 
and Manohar's RPI, to order 1/m?. 

In order to understand the nature of the difficulty and how to resolve 
it, let us briefly review the essentials of the HQET procedure. A priori the 
correct procedure for calculating an amplitude to some order in 1/m (and in 
as) is to compute the relevant Feynman diagrams of full QCD and to expand 
the result in powers of 1/m. HQET works essentially in the reverse order, 
which is the key to its conceptual and computational advantages. Having 
separated out the fixed large momentum as in eq. (1), the QCD vertices 
and propagators are first expanded in powers of 1/m and then used to 
compute Feynman diagrams. Because of the different ultraviolet divergence 
structures that arise depending on whether one expands in 1/m first or last, 
this naive version of HQET only agrees at tree-level with full QCD. However, 
the two procedures can be made to agree ( "match" ) by adding local terms to 
the HQET lagrangian. Once the extra terms in Chqet are determined (to 
any particular order) by matching to some number of full QCD amplitudes, 
they are universal, and can be used to compute other amplitudes without 
further reference to full QCD. This procedure is quite general and applies 
even if the quarks were spinless. In the spinless case there is no difficulty 
in tracking w-dependence, deriving the form of RPI for the effective theory, 
and enforcing the constraint on the form of the effective lagrangian |^ . 

For spin-1/2 quarks there is an extra complication. When the quark 
propagator is expanded for <^ m m eq. (1), one finds 

' ' ^Oih (3) 



— m + ie 2 k.v + ie m' 

This shows that in the heavy quark limit m — > oo, not all four spinor degrees 
of freedom propagate, but only the two components projected out by (1 -|- 
f) /2. The usual practice is to interpolate these spinor modes in HQET with 
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a heavy quark field i^+v, 

^V+i) = ip+v (4) 

(For an interpolating field this makes sense even when m is large but not 
infinite.) Of course the expansion of eq. (3) is invalid in loops, where k need 
not be small, and all four spinor degrees of freedom are important. Let us 
denote the remaining degrees of freedom by ip-v, 

fip-^ = -ip-v (5) 

We can refer to if^j^^ and V-ij as the "particle" and "anti-particle" fields 
respectively, since eqs. (4,5) are just the Dirac equations for positive-energy 
particles and anti-particles with fixed momentum mv. 

We can understand the status of V'-i) by observing that for A;^ <C m in 

(1), 

2 i)-m + ie 2 2 2m ^m?'' ^ ' 

We see that relative to the ground state of the heavy quark sector, ip^^ is 
a massive mode, with mass gap 2m. Like any heavy mode in field theory 
we can integrate out its virtual effects and omit the field from the effective 
theory. This is what is done in standard HQET, so Chqet depends only 
on and the light quarks and gluons. Normally in field theory heavy 
and light particles transform independently under symmetries. Integrating 
out the heavy particles does not alter the symmetry properties of the light 
particles. RPI is different, in that under eq. (2), 'ip±^ clearly mix. Integrating 
out the heavy mode tp-^ necessarily complicates RPI in the effective theory. 

In this paper, I derive RPI from first principles, in two versions of HQET. 
In section 2, 1 describe a version of HQET in which the anti-particle field ip—v 
is explicitly present, but as a non-propagating auxiliary field.[| As in stan- 
dard HQET, the 1/m expansion is performed before computing diagrams. 
The fact that both are present makes RPI as simple as it is in the case 
of heavy spin-0 particles. One can in fact write the most general form of the 
RPI-constrained effective lagrangian. This is done in section 3. In the sec- 
ond step of the program, described in section 4, the ip-^ field is completely 
integrated out, resulting in the general form of the RPI-constrained effective 
lagrangian in the standard HQET formulation. In section 5, I check some 
of the consequences of RPI on Chqet at the leading orders of 1/m. Section 
6 contains some discussion and comparison with earlier proposals. 

t ■!/)_„ plays a role in simplifying RPI analogous to the role of auxiliary fields in super- 
symmetric field theory. 
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For simplicity, in this paper I only consider HQET applied to the case 
of a single heavy quark in interaction with light quarks and gluons. While 
the character of the arguments and results presented in this paper are non- 
perturbative in nature, the proofs are given within perturbation theory, in 
the context of minimal subtraction and dimensional regularization. I will 
address the topic of RPI in non-perturbative lattice HQET elsewhere. 



2 HQET with an Auxiliary Field 

To focus on the residual momenta as in eq. (1), we make the standard 
field redefinition, 

^4^)^e™V(^), (7) 

where ip is the full QCD quark field and ip^ will be our HQET field, and m 
is the heavy quark mass. We can further decompose ip^, 

i^±v = -Y^^v. (8) 

There is an ambiguity in just what is meant by the quark "mass". The 
canonical choice is to choose m so that the residual mass term for ipj^^ van- 
ishes in the HQET. In general though, any choice is allowed if it makes the 
residual mass ^ m (see ref. ^J). The minimal subtraction mass parameter 
of full QCD will be separately denoted by m(/i). 

At tree- level, expanding in 1/m before or after computing Feynman dia- 
grams makes no difference, so matching the HQET to full QCD only amounts 
to accounting for eq. (7), 

= 'ijj_^^{iD.v + m — m)'>p+v 
— 'ip-yijn + fh + iD.v)ip-v 
+ i^+vip±^-v + Tp.JpxiJ+v, (9) 

where ]p± = p — D.v^. The gauge field and light quark terms have been sup- 
pressed. The canonical choice for m at this order is clearly m = fh. In using 
this lagrangian the HQET approach tells us to expand in 1/m (1/m) before 
computing diagrams, so we see that the V'-d propagator in a background 
gauge field is really completely local, 

1 1 iD.v , , 

+ 7Z-^-- (10) 



m + m + iD.v m + rh {m + rh)'^ 
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Order by order in 1/m this auxiliary field ip-^ can be integrated out to yield 
interactions for the propagating field. 

At the quantum level, matching corrections to Chqet are non-trivial 
because in full QCD we expand in 1/m after computing diagrams while in 
the HQET we expand the Feynman vertices and propagators in 1/m before 
computing diagrams. The two steps do not commute inside loops because of 
ultraviolet divergences, but since such divergences are local, order by order 
in the loop expansion we can compensate by adding purely local terms to 
the HQET lagrangian. To state the matching conditions precisely, we add 
sources to full QCD and to the HQET, 

SCqcd = 7?e^'""-"V + h.c. 
SChqet = rjV'i, + h.c, (11) 

and demand that the corresponding Green functions agree order by order 
in 1/m and a^.^ 

The most general form for the HQET lagrangian is, 

Chqet = + f + J] TW.a^")]^'. (12) 

n>2 

n>2 

+ ^_J5-t;.F+5]T(").cT(")]^_, 

n>2 

n>2 
n>2 

where S, V^, T^\,^„ are local hermitian operators constructed from the light 
fields, covariant derivatives, and v, and where cj(") is the totally antisym- 

■'■These could be either connected or IPI Green functions. Note that normally the 
relation between IPX vertices in full and effective theories is complicated when a massive 
mode is completely integrated out. The reason is that a diagram in the full theory that 
can be cut in two by cutting an internal massive mode propagator can become IPX in the 
effective theory where the massive mode is absent. In the present discussed in the 

introduction, the massive mode is i/>_„ , but because we have left it in the effective theory 
as an auxiliary field, we have a simple equality of the IPX effective actions of the full and 
effective theories after matching. 
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metrized product of n (4 + e)-dimensional 7-matrices, 



7[/.i7/^2"-7/.„i; n = 4/c - 1,4/c - 2, 

7[/.i7/^2"-7m„]; n = 4k,Ak-3. (14) 

In four dimensions u*^'^^ and fx^^^ can be reduced by introducing 75, but I have 
refrained from doing this to avoid problems of dimensionally continuing 75. 
In four dimensions we also have a^"^ = for n > 5, so they correspond to 
evanescent operators, only occurring in the counterterm lagrangian of the 
effective theory.^] 

The absence of operators higher than bilinear in the heavy quark field 
is because such vertices make no contribution to (dimensionally regulated) 
amplitudes in the effective theory describing a single heavy quark. 

The different Dirac structures are of different orders in 1/m, which can 
be worked out by finding the lowest dimension operators that can occur in 
HQET with these structures, and respecting QCD symmetries. However 
one of these symmetries is RPI itself, which is discussed next. This power- 
counting exercise is therefore deferred till the end of the next section. 



a 



(n) 

Hi. ..fin 

(n) 

Hi. ..fin 



3 Reparameterization Invariance 

The HQET with the auxiliary anti-particle field satisfies a very simple form 
of RPI, Chqet being invariant under 

v^v + 5v- i)^{x) e'"'^''-='ip^{x). (15) 

At tree-level for example, this invariance is manifest in the first line of eq. 
(9). Beyond tree-level the proof of RPI goes as follows. By eq. (7), the quark 
field of full QCD, iIj{x), is invariant under eq. (15). This means Cqcd is 
invariant except for the source term, eq. (11). Clearly the QCD partition 
functional then satisfies, 

^The reason the evanescent operators only appear among the minimal subtraction 
counterterms is just the generalization of the arguments of ref. Q . It is possible that the 
technical problem of evanescent operators can be avoided by using dimensional reduction 
as the regulator. See for example the discussion in ref. [^. I have not investigated this 
possibility in the present context. 
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where is the gauge field source. This relation holds even when the func- 
tional integral is dimensionally regulated and minimal subtraction countert- 
erms added. (Dimcnsonal regularization is a "good" regulator for reparam- 
eterizations, cq. (15).) After matching, this RPI and renormalizcd partition 
functional is equated (to whatever order in ag and 1/m one is working) to 
the renormalized HQET partition functional, 

Z[ri, fj, Ji,, V]= J VA^Vi^^V^J^'^'^^'^ [CHQET+C,.t.+J.A+r]^.+ii,ri\ ^ 

where Cc.t. denotes the HQET minimal subtraction counterterms, and the 
Fadeev-Popov ghost determinant is subsumed into the gauge field measure. 

Noting that the source term, eq. (11), and the regulated measure are both 
invariant under r/ e^"^^^'^r],tjjjj — >■ e^"^^^'^tpy, we see that 

VAnVljjyThlj I ICHQET+SCHQET+^c.t.+SCc.t.+J-A+rftjJv+i'yV] ^(^Ig^ 



where, 

SC = C{e'^'---^,,-^,e-'^'---,A^, v + 5v) - C{il^,Av. v). (19) 
Eq. (16) therefore reads, 

j PA^PV^P^^ j dx{5CHQET+5Cc.t)e'i''^'' CHQET+c,.t.+J.A+v^^+^,v = 0. 

(20) 

It follows that 

J dx{6CHQET + SCc.t.) =0, (21) 

since all insertions of this operator into Green functions vanish. In minimal 
subtraction we separately have, SChqet = and SCc.t. = since they have 
distinct 1/e-dependence. This proves the RPI of Crqet- 

The RPI constraint on Chqet has a simple solution: and can 
only occur in the lagrangian with ipy,ipv (see eq. (12)), in the combina- 
tion + iD^/m. (This is only true for covariant derivatives acting on 
the heavy fields, not on the light fields, which are trivially RPI.) That is, 
S", V^j,T(") arc constructed out of + iD^/m. The proof is straightfor- 
ward. Instead of working with ipy,ipv directly, we can write any effective 
lagrangian in terms of the combinations tp^e^^'"'^ and e~^"^'"'^'ijjy which, as 
we noted earlier, are invariants under eq. (15). A general term of Chqet 
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has the form ip.^e^^^'^ f[v^^Dy)e~^^^'^ipv where / is some monomial in its 
arguments (and can be a matrix in spinor, Lorentz or color spaces), mul- 
tiplied by a term involving only light fields and their covariant derivatives. 
In this decomposition only the argument of / transforms under eq. 
(15). Furthermore the only invariant combination of u^'s alone is trivial, 
= 1. Therefore imposing RPI just translates into the ^-independence of 
/. The phase factors can be cancelled against each other by noting that 
Dfj^e^"^^-^ ... = e^"^^'^{imv^ + D^)... . Thus, the general RPI term in Chqet 
has the form il)^f{imv^ + D^)ip^, multiplied by a factor consisting of light 
fields, thereby completing the proof. We see that in the auxiliary field for- 
mulation, RPI is as simple as in theories with spin-0 heavy quarks Q. As 
emphasized by Luke and Manohar, the central reason for the complication 
of RPI in the standard formulation of HQET is that the heavy quark field 
must satisfy a f-dependent constraint. By including the auxiliary field in 
the present formulation, there is no constraint on ip^j = ipj^^ + ip-v- 

We are now in a position to give the power-counting for the operators 
of the various Dirac structures, described at the end of the last section. 
They are restricted by conventional QCD symmetries, and by RPI. The 
must have at least one pair of anti-symmetrized Lorentz indices to 
be contracted with the cr^"'^ These indices can only arise from a pair of 
reparameterization-covariant derivatives, Vp^+iD^/m. That is, T*^") = ...{v+ 
iD/m)[p...{v + iD/m)^j... . However, terms in which any u's are pulled out 
of the parentheses do not survive the anti-symmetrization unless there are 
also derivatives between the two sets of parentheses. It follows that the 
r(") operators must contain at least two derivatives, D^^. So by dimensional 
analysis the T^") can be at most 0{l/m). The lowest dimension operator 
that can appear in S is the unit operator, with a coefficient of order m 
by dimensional analysis. Similarly the lowest dimension operator that can 
appear in is t;^ + iDf^/m, again with a coefficient of order m. There is 
just one subtlety which can already be seen at tree-level. Though S and 
Vp are both of order m, S + v.V ~ 0{1) if we choose the residual mass to 
vanish (or be at most 0{Aqcd))- To summarize, 

S -v.V 0{m); S + v.V^O{l); fx ~ 0(1); T^") < 0(l/m). (22) 

4 Integrating out the Anti-Particle Field 

The auxiliary field formulation of HQET is a valid calculational scheme 
in its own right, once matched to full QCD. However, having solved the 
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RPI constraint in this formulation we are free to integrate out the auxihary 
field and return to the standard formulation of HQET in terms of 
(once we turn off the source for tp-v)- The result will be the general RPI- 
constrained form of the standard HQET lagrangian. Since ip-^ occurs at 
most quadratically in eq. (13) we can easily integrate it out to get, 

n>2 n>2 

X [S- v.V + ^^"^•^^"V--)!^^ + r^"^-^^"^]}V'+.. (23) 

n>2 n>2 

This expression requires some explanation. The (—v) subscript appearing 
in the operator inverse is an instruction to perform the inversion completely 
within the (1 — ^)/2 subspace, not the full spinor space, thereby getting a 
result that also lives in the subspace. Secondly, this inverse is really local. 
The reason is that the dominant part oi S — v.V is a constant of order m. 
Therefore expanding the operator inverse in powers of 1/m yields a series of 
local operators, similarly to eq. (10). The last observation concerns renor- 
malization. When integrating ip^^ out of the auxiliary field formulation, the 
effective lagrangian in the functional integral contains counterterms. There- 
fore each operator Oq = S, V, T^") appearing in eq. (23) has the form, 

Oo = Oifi)+ Opoles, (24) 

where Opoies is a sum of l/e" pole counterterms obtained in the auxiliary 
field formulation of HQET. The standard HQET lagrangian is just the finite 
part of eq. (23), the poles providing the counterterms for the standard 
HQET formulation. Clearly, given eq. (24), to all orders in ag and 1/m this 
finite part must also have the form of eq. (24), where the operators are the 
0(/i).p] That is, both the renormalized and unrenormalized standard HQET 
lagrangians have the form of eq. (23). RPI restricts the 0{fi) and the Opoies 
to be constructed from + iD^/m. This result agrees with ref. @. 

We can view the form of eq. (23) as the solution to demanding invari- 
ance under a reparameterization transformation of ipj^v ■ We can derive this 

^The products of (4 + e)-dimensional 7-matrices encountered in eq. (23) can be re- 
expressed as linear combinations of the fully anti-symmetrized products, l,7^,cr'"\ and 
simplified using (l-|-^)/2-projected identities to arrive at a canonical form for the standard 
HQET lagrangian. These particular manipulations do not introduce any explicit powers 
of e, which would have changed the separation of the effective lagrangian into the finite 
part and minimal subtraction counterterms. 
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transformation by noting that the Gaussian integral performed to ehminate 
is the same as using the ijj-^ "equation of motion" following from eq. 

(13), 

= -[S-v.V+J2 r("^c7(")](-i^) [Y± + E (25) 

n>2 n>2 

(See the remarks after eq. (23) regarding the interpretation of the operator 
inverse). Thus the reparameterization transformation in the auxiliary 
field formulation corresponds to the following transformation in standard 
HQET, 

W+v' — — ^ — e Wv 
~ 2 ^ 

n>2 n>2 

r-, • e 

= |1 + imov.x + ^ — y. 

[s -v.v+Y. r(").a(")]-i^) [n + E (26) 

n>2 n>2 

where 0{5v'^) terms have been dropped in the last line. This general form 
also agrees with Kilian and Ohl @. However, the specific forms of 5, V, T^") 
must be determined by matching the auxiliary field formulation of HQET 
to full QCD. It is straightforward to see that the operators appearing in eq. 
(26) are either the 0(/x) or Oq depending on whether one is considering RPI 
of Chqet or Chqet + Cc.t. respectively. 



5 Examples in Standard HQET 

It is a well-believed expectation that RPI should constrain the term in the 
standard HQET lagrangian, — 2^'i/;_,_^D^V'+t;) to have unit coefficient to all 
orders in (for example the analogous statement is straightforward to prove 
for heavy scalars ) . Let us see how this emerges from the general form, eq. 
(23). One can easily check that the T*^") can have no effect on the 0(l/m) 
operator of interest. Thus the part of the lagrangian with the right Dirac 
structure to this order is given by 

Chqet = ^^+,{8 + v.V -Y^S - v.V)-^Yi_}i>+, + ... , (27) 
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where S" is a function f{v + iD/m) )) = /(I + 2iv.D/m — D /m ), and 
= g{l + 2iv.D/m — D"^ /m?){v^ + iZ)^/m). Actually, RPI permits more 
general forms for S and V which can be reduced to the above by using the 
fact that + iD^/m commutes with + iD^/m up to terms involving the 
gauge field strength, such terms being discarded into the ellipsis above. To 
order 1/m the general forms of S and are, 

D2 (v.DY 

S = -mi + Ziiiv.D ) + ki- — 

2m m 
n2 (y r)\2 Tj 

= [m2 + Z2{iv.D-—) + k2^-^]{v^ + i^). (28) 
Zm m m 

Substituting into eq. (27), 

C-HQET = V'+„{"T'2 - mi 

+ Zi + Z2){w.D -—) + —iv.D - ^ + ... 29 

2m m m^(mi + m2j 

It remains to put Chqet into canonical form. The m2 — mi term is the 
"residual mass" term. We choose m to make the residual mass vanish, 
mi = 771-2. Next we perform wavefunction renormalization to ensure that 
ipj^^iv.Dtpj^y has canonical unit coefficient. The result is 

Chqet = ip+^iiv.D - ^)i^+v + ••• , (30) 

as we wished to prove. 

A second example involves spin-dependent terms in the standard effective 
lagrangian. Both Chen's and Luke and Manohar's versions of RPI predict 
that the coefficients of the operators, 

O2 = ^?+,c7"^7;'^{I?„,G^.}V+., (31) 

are related by 

'2Cmag = 6*2 + 1 (32) 

(when the heavy quark kinetic term and residual mass terms are put into 
canonical form) 0]. Ref. 0] showed this to be true independently 
of the earlier RPI proposals. It is therefore of interest to know whether 
eq. (32) follows from our general form, eq. (23). The answer is yes, as 
a straightforward but tedious calculation (along similar lines to the last 
example) shows. 
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6 Discussion 



At tree-level, S = —ifi, = mv^ + iD^, T^") = 0, and if we choose m = m, 
eq. (26) becomes the transformation proposed by Chen 0. Beyond tree- 
level, there are definitely matching corrections to 5, V, T^"^ in the auxil- 
iary field formulation of HQET, so the reparameterization transformation of 
standard HQET, eq. (26), is necessarily corrected beyond Chen's tree-level 
form. This means that Chen's proposed form of the RPI transformation 
must be wrong. Therefore at first sight, it is surprising that Finkemeier, 
Georgi and Mclrvin Q could prove all the predictions of Chen's form of 
RPI to 0{l/m?) in Chqet, such as eq. (32) for example. We are in a posi- 
tion to understand this. The results to this order in the effective lagrangian 
are completely determined by the form of the reparameterization transfor- 
mation to order 1/m, so that the T^") are irrelevant inside eq. (26), by eq. 
(22), and we can use the forms in eq. (28) for S and V. When we enforce 
the vanishing of the residual mass, mi = m2, we find, 

V'+^z = {1 + im6v.x + ^ + ^ ^}^+. + 0{l/m^). (33) 

2 2 2m 

This is precisely the same as Chen's transformation to order 1/m, thereby 
explaining the success of Chen's proposal to order 1/m^ in the effective 
lagrangian. Beyond this order, eq. (26) and Chen's transformation no longer 
coincide, and only eq. (26) is correct. 

Eq. (26) is also not the RPI transformation proposed by Luke and 
Manohar . However to order 1 /m^ , ref . has shown that there is a field 
redefinition (which does not affect the S-matrix) which, when compounded 
with the Luke-Manohar transformation, yields eq. (33). It is possible this 
agreement with eq. (26) up to field redefinitions persists at higher orders in 
HQET. See ref. |^ for some more discussion. 

The auxiliary field formulation of HQET is a viable calculational scheme 
which makes RPI extremely simple, given by eq. (15). If one prefers the 
standard formulation of HQET, the correct form of RPI is the one obtained 
by Kilian and Ohl Q and rederived in this paper. 
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